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L'invention est relative k un proc€d6 pour determiner localement la forme des 
horizons geologiques k partir des donn^es obtenues pai- traitement de mesiires 
sisnaiques. 

L'invention est €galement relative k un dispositif pour determiner localement 
la forme des horizons gdologiques k partir de donn€es obtenues par traitement de 
mesures sismiques. 

L'invention est enfin relative k un produit-programme d'ordinateur permettant 
le fonctionnement d'un dispositif programmable selon l'invention, pour la mise en 
oeuvre d'un proc^dd selon l'invention. 

Le document US 6 151 555 est relatif k un proc6d6 de traitement de signaux 
sismiques et ^ un dispositif pour la determination d'un cube de valeur de variance. Ce 
procede prevoit de calculer pour chaque noeud d'une matrice sismique 
tridimensionnelle une valeur de variance dgale a la variance des amplitudes 
sismiques au voisinage du noeud considere. Le voisinage du nceud consider^ est un 
sous-ensemble des noeuds de la matrice tridimensionnelle dont la distance au noeud 
consider^ est inferieure a un seuil choisi. Ce precede repose sur I'hypoth^ise que la 
valeur de variance ainsi definie est importante pour les noeuds proches d'une 
discontinuite. 

Le document WO 99/64896 decrit un procede d'inteipretation de donnees 
sismiques, particuliSrement adapte pour identifier les limites de reservoirs 
d'hydrocarbures. 

Un premier but de l'invention est d'ameiiorer les analyses du profil des 
horizons geometriques. 

Un deuxifeme but de l'invention est de feciliter la determination des vectevirs 
unitaires normaux aux horizons d'un cube sismique tridimensionnel. 

Un troisi&me but de l'invention est de faciliter la detei-mination d'un indice de 
courbure des horizons dans un cube sismique tridimensionnel. 

L'invention a pour objet un procede pour determiner localement la forme des 
horizons geologiques, dans lequel on construit iine fonction continue SjjjcCt) par 
interpolation ou approximation des traces sismiques discretes d'une matrice sismique 



tridimensionnelle, ladite fonction etant d€sign6e comme "trace sismique locale 
continue", comportant les Stapes suivantes : 

a) r nmmft d^A^age optimal de deux traces sis miques locales 

continues voisines, la valeur de ddcalage rendant maximale leur fonction de 
5 correlation ; 

b). retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centrale*' Sy,k(t) le sous-voisinage consistant en des traces voisines Spq,k(t) 
coiTespondant k des decalages optimaux associds i des correlations Rij,pq,k(h) 
supdrieures k un seuil predetermine compris entre 0 et 1 ; 

10 c). definir pour chaque trace sismique locale continue Spq^(t) du voisinage 

conditionnel une valeur de lesidu par rapport k ladite trace sismique locale continue 
"centrale" Sij,k(t) comportant des coefficients parametriques ; 

d). determiner les coefficients parametriques par nainimisation d'ensemble des 
residus sur le voisinage conditionnel. 

15 Selon d'autres caracteristiques de Tinvention : 

- les valeurs de residu de retape c) sont des valeurs absolues de forme 
differentielle parametrique, par exemple Aj,pq.k = I ^iijn .(p-i) + bij.k -(q-j) - hij.pq;c I » 

oil aij,k et bij^ sont les paramfetres et hij,pq^ est le decalage optimal retenu entre 
la trace sismique locale continue Sij,k(t) "centrale" et une trace sismique locale 
20 continue Spq,k(t) voisine, 

- la minimisation d'ensemble des residus p,^^^^ sur le voisinage conditionnel 
comporte une minimisation d'une sonome de puissances des residus du type suivant 

C ^ (i,j,k) = X(^u.pq.k^" » puissance superieure ou egale k 1, 

Le cas echeant, chaque terme (pij.pq.k )"de cette somme pent Stre ponderee par 
25 un facteur fonction croissante du coefficient de correlation flj,pq.k(hij.pq,k)- 

- les coefficients ay^k et byj^ determines k Vetape d) sont utilises pour definir un 
vecteur unitaire N(i j,k) de coordonnees N\iJ,k) , N^(ij,k) , N'(i,j,k) , par exemple 
de la forme: 



N''(id.k) = -= 

Ny(ij,k) = , ^'^'^ = 

N*(ij.k) = 



VCaij,k)^+(bij,k)^-l 
et on choisit ce vecteur unitaiie comme vecteur unitaire «nonnal» orthogonal 
k I'horizon au nceud (i, j , k), 

- on utilise I'ensemble des vecteurs unitaires nonnaux N(i, j, k) pour d^finir le 
champ des vecteuts unitaires nonnaux, en vue de I'afFichage sur dcran des profils des 
horizons sismiques, 

- on choisit comme indice de courbure une fonction croissante des r6sidus 
comportant les coefficients param^triques d€tennin€s h I'dtape d). 

L'invention a dgalement pour objet un dispositif pour la mise en ceuvre du 
proc€d^ selon l'invention, comportant des moyens pour utiliser comme d^calage 
optimal de deux traces sismiques locales continues voisines la valeur de ddcalage 
rendant maximale leur fonction de corr^ation, des moyens pour retenir comme 
voisinage conditionnel d'une trace sismique locale continue "centrale" Sijjt(t) le sous- 
voisinage consistant en des traces voisines Spq,k(t) correspondant k des d&alages 
optimaux associes k des correlations Rij,pqjc(h) sup^rieuies k un seuil pr^ddtermme 
compris entre 0 et 1, des moyens pour d^finir pour chaque trace sismique locale 
continue Spqjc(t) du voisinage conditionnel une valeur de r&idu par rapport k ladite 
trace sismique locale continue "centrale" Sijj,(t) comportant des coefficients 
param^triques et des moyens pour determiner les coefficients param€triques par 
minimisation d'ensemble des r^sidus sur le voisinage conditionnel. 

Selon d'autres caract6ristiques de l'invention : 

-le dispositif comporte des moyens de memorisation et des moyens de 
visuaUsation de parametres sismiques determines k I'aide du precede selon 
l'invention. 



L'invention a enfin pour objet un produit-programtne d*ordinateur, 
comportant des elements de code de programme pour executer les etapes d'un 
procddd selon Tinvention, lorsque ledit programme est ex^cutd par un ordinateur. 

L'invention sera mieux comprise grSce k la description qui va suivre donnfe k 
5 titre d'exemple non limitatif en reference aux dessins annexes dans lesquels : 

- La figure 1 reprdsente schematiquement une matrice sismique 
tridimensionnelle et une trace sismique locale continue. 

- La figure 2 reprdsente schematiquement, un exemple de voisinage local Nij,k 
consistant en un ensemble de trace sismique locale continues Spq^(t) voisines de la 

10 trace centrale de r6fi^rence Sij.k(t) elle mSme situ€e a la verticale du g^ophone Gy et 
interpolant les amplitudes sismiques autour de t=k. 

- La figure 3 repr&ente schematiquement, une section verticale de voisinage 
local conditionnel Nij.k(r) d'une trace sismique locale continue de rdKrence Sij.k(t)» 

- La figure 4 repr&ente schematiquement, une section verticale dhine matrice 
15 sismique tridimensioimelle avec un ddcalage optimal hij,pqjc et un voisinage 

conditionnel Nij,k(r). 

- La figure 5 repr&ente schematiquement, une section verticale d'une matrice 
sismique tridimensionnelle analogue k la figure 4 avec presence d'une faille ou 
discontinuity. 

20 - La figure 6 reprfeente schematiquement, un schema fonctionnel d'un 

procede selon l'invention. 

En reference aux figures 1^5, les elements identiques ou fonctionnellement 
equivalents sont designes ou reperes de manifere identique, 

Sur la figure 1, une matrice sismique tridimensionnelle est obtenue par releve 
25 de mesures enregistrees par des geophones Gy disposes sur un reseau x,y en des 
points de coordonnees ij. Le releve des mesures echantillonnees dans le temps est 
reprdsente suivant un axe t descendant representatif de la profondeur ou d'une 
verticale descendante a partir de la surface de la terre ou de la mer. Les mesures sont 
caracterisees par leur amplitude, par exemple une amplitude relevee par le geophone 



G'J au temps ou k la profondeur d'dchantillonnage tk. La mesure discrke effectu€e par 
le g&phone Gy au temps ou k la profondeur tk est appelfe amplitude sismique Sij,k. 

L'ensemble des amplitudes sismiques correspondant h. un g&>phone Gy de 

coordonnfes ij est une matrice unidimensionnelle (Syi, Sij2. Syk,...., Syw) 

appeI6es trace sismique discrete, car cette matrice unidimensionneUe correspond k la 
trace selon le point de coordonn&s horizontales iJ de la matrice sismique 
tridimensionnelle obtenue par mesures sismiques. 

L'axe vertical t orient^ selon une verticale descendante d€signe usuellement le 
temps, mais peut dgalement 6tre traits pour repr&ent6 une profondeur k partir de la 
surface. 

L'invention conceme aussi bien I'application k une troisifeme coordonn^e t 
representative du temps, qu'k une troisifeme coordonnde t repr&entative de la 
profondeur. 

A partir de la trace sisnaique discrete sitafe k la verticale d'un ^ophone Gy 
on definit, par intetpolation ou approximation des valeiu^ discretes autour de t?=tk=k, 
une fonction continue S|j,k(t) qui est d^signfe comme "trace sismique locale 
continue". Les mdthodes d'approximation ou d'interpolation de valeurs discretes pour 
engendrer une fonction continue sont nombreuses, et comprennent notamment les 
interpolations ou approximations polynomiales, ainsi que les interpolations ou 
approximations trigonom^triques polynomiales. 

Toute autre vaiiante d'interpolation ou d'approximation foumissant une 
fonction continue peut egalement Stre appliqu€e k la px^ente invention pour fourair 
une "trace sismique locale continue". 

Sur la figure 2, plusieurs traces sismiques locales continues d^finissent un 
exemple de voisinage d'une trace sismique locale continue "centrale" Sij,k(t) de 
reference. Le voisinage d'une trace sismique locale continue Sij,k(t) est defini comme 
l'ensemble des traces sismiques locales continues dont les indices spaciaux 
horizontaux p,q sont voisins des indices spaciaux horizontaux iJ de la trace sismique 
locale continue de r6f(5rence. 

A titre d'exemple, les coordonndes spatiales horizontales p,q correspondant au 
g^ophone Gpq sont voisines des coordonnfes spatiales horizontales iJ correspondant 



au g^ophone Gy si les valeurs absolues des differences i-p et j-q sont infdrieures k des 
entiers donnas, par exemple k 2, 

DaasxacaSyjCQmmejgprfeentd^^a figure 2, la trace sismique loralft confimie 

Sijjc(t) est associee h huit traces sismiques locales continues voisines entoui^ant la 
5 trace sismique locale continue "centrale" Sij,k(t). 

Dans le cas de traces sismiques locales continues produites par des mesures 
sismiques, la forme des horizons geologiques introduit des ddcalages verticaux entre 
des traces sismiques locales continues voisines. Dans le but de ddteroiiner des 
relations entre deux traces sismiques continues locales voisines Sij,k(t) et Spq,k(t) 
10 centr^es sur la m§me coordonnde verticale d*€chantillonnage t=tk=k et correspondant 
& des coordoxmdes spatiales i J et p,q differentes, on calcule la fonction de correlation 
Rij,pq,k(h) des deux traces sismiques locales continues. 

La fonction de correlation de deux traces sismiques locales continues 
voisines est obtenue par la formule suivante 



15 %pq,k(h) = 



Cij,P.,ic (h) 



VCu.u..(0)-Cp..p,..(0) 
ot le numerateur correspond k la fonction de covariance de SyjcCt) et Spq,k(t) 

Ifc+A 

obtenue par I'expression suivante Cij,pq,k(h) = J Sy j, (t) . S (t + h) . dt . 

tk-A 

Dans cette int^grale definissant Cij,pq,k(h) , le param^tre A d^finit une «fen8tre 
20 verticale d'investigation » autour de t=tk=k. 

Par exemple, si les traces sismiques locales continues SyjcCt) et Spq,k(t) sont 
des polynSmes trigonometiiques de la forme suivante interpolant les donn6s 
sismiques 



Siik(t) = 5^af^.cos(sa)t) + b3J**'.sin(sfl)t) 

m 

25 Sij.k(t) = 51aij''.cos(sfl)t) + b[^*.sin(sa)t) 



alors on peut montrer que la fonction de covariance Cij,pqjc(h) est elle meme un 
polynCme trigonom^trique de la fonne suivante oil les coefficients As'^^^^'^ et Bs'-**^*^ 
dependent des coefficients as'^*^, a/^\ hs^*^ et bs^*^ : 



Un rdsultat math^matique classique et connu est que la fonction de correlation 



lorsque cette fonction de correlation se rapproche de 1. 

L'^tude des fonctions de correlation des traces sismiques locales continues 
Sij,k(t) et Spq,k(t) permet de d^finir conime ddcalage optimal la valeur hy.pq^k de h 
correspondant au maximuna de la fonction de correlation Rij,pq,k(h), c*est-a-dire au 
maximum de la fonction de correlation le plus proche de 1. Contrairement a Tart 
existant base sur une formulation discrete de la fonction de correlation Rij,pq,k(h), 
Tutilisation d'une formulation continue des traces sismiques et done de la fonction 
de correlation permet d'obtenir un decalage optimum qui n'est pas contraint k Stre un 
multiple entier du pas d*echantillonnage suivant Taxe vertical correspondant k la 
variable t , 

La definition du decalage optimal hy^pq^k de deux traces sismiques locales 
continues voisines permet d'obtenir une premiere approximation des horizons passant 
par cette trace sismique locale continue de reference Sy^kCt), comme represente aux 
figures 3 a 5. Contrairement a Tart existant, le fait d'utiliser des decalages optimaux 
qui ne sont pas contraint a 8tre un multiple entier du pas d*echantillonnage suivant 
Taxe vertical permet d*eviter des erreurs connues sous le nom de «aliasing». 

Sur ces figures 3 ^ 5, on definit un voisinage conditionnel Nij.k(r) de la trace 
sismique continue locale de reference Sy^kCt) qui est un sous-voisinage du voisinage 
initial des traces sismiques locales continues retenues pour le calcul de correlation et 
de definition des decalages optimaux. Le voisinage conditionnel Nij,k(r) est choisi de 
maniere que, pour toute trace sismique continue locale Spq,k(t) appartenant h Ny^kCr), 
le decalage optimal hij,pq.k correspondant est tel que la correlation Rij,pq,k(hij.pq.k) entre 
Sij,k(t) et Spq,k(t) est superieure k un seuil predetermine r compris entre 0 et L 



m 

Cij.pqjc(h) = ^ AfP^^cosCsQ) h) + B;j•P^•^sin(sa) h) 




<^ij.pq.k (h) 



traduit une similarity des fonctions con^l^es 



Comme on le voit sur la figure 5, la fixation du seuil r permet egalement 
dfengiober dans un voisinage conditionnel une discontinuite ou une faille g^ologique, 
ce qui constitue un avantage important par rapport h Tart antdrieur, 

Les dispositions precitees foumissent ainsi une modelisation continue 
5 permettant la mise en oeuvre de Tinvention. 

Sur la figure 6, la modelisation numerique d^crite en reference aux figures 1 k 
5 est effectute k T^tape 1 00 pour mettre en ceuvre un proc€de selon Tinvention. 

Apr^ cette 6tape 100 de modelisation num€rique, on d^finit k Tdtape 101, 
pour chaque trace sismique locale continue Spq,k(t) du voisinage conditionnel NyjcW 
10 une valeur de r&idu par rapport k ladite trace sismique locale continue Sij,k(t) 
centrale de rdfdrence comportant des coefficients paramdtriques. 

La valeur de rdsidu Pij.pqJc choisie traduit Tecart entre Thorizon k la verticale 
du geophone Gpq et le plan tangent a Thorizon a la verticale du gdophone Gy. 

A titre d'exemple non limitatif, on pent d^finir comme valeur de r^sidu une 
15 valeur absolue de forme differentielle paramdtrique, du genre : 

pij,pq.k = laij,k.(p-i) + bij,k.(q-j) - hij,pq.k|, oii ay^ic et bij,k sont des paramfetres et 
hij,pq,ic est le ddcalage optimal retenu entre la trace sismique locale continue Sij,ic(t) 
"centrale" et une trace locale sismique locale continue Spq,k(t) voisine. 

Aprfes avoir d^fini Tensemble des residus sur le voisinage conditionnel Nij,k(r), 
20 on determine les paramfetres ay^k ct by^ qui correspondent k une mininnisation 
d*ensemble des residus Py^pq^k sur le voisinage conditionnel k Tetape 102. 

Le critfere de minimisation choisi est arbitraire, et pent par exemple comporter une 
minimisation d'une sonmie de puissance sup6rieure k 1 des residus, de la forme 

C"(i J ,k) = ^ (Py.pq,k T , ou n est une puissance r€elle positive superieure k 1 , 
25 Le cas echeant, chaque terme (Pi|,pq,k )"de cette somme pent dtie pond^ree par 

un facteur fonction croissante du coefficient de correlation jF5j.pq.k(hij,pq.k)- 

La minimisation effectuee k Tetape 102 produit comme r€sultat a Tdtape 103 
la meilleure estimation des coefficients ^j,k et by^c qui sont les coefficients directeurs 



9 



de rhorizon passant par le point de coordonn^es horizontales ij et de coordonnde 
temporelle ou de profondevir k. 

A r^tape 104, selon una premiere variante de I'invention, on d€tenmne les 
vecteurs unitaires normaux & I'horizon passant par le point de coordonn&s spatiales 
i, j et de coordonnde temporelle ou de profondeur k par des formules de type 
connu : 

N''(ijjc)= 



VKJ' + (b,.,)^ + l 



Ny(ij,k)= ^'j" 



N'(iJ3c) = 



1 

>/(aij.k)^ + (bij.k)^-l 



Le vecteur unitaire normal ainsi obtenu ou tout vecteur qui lui est parallele 
permet une visuaUsation sur ^cran plus efficace que les visualisations de I'art 
ant^rieur, de mani^re k reveler les paramfetres sismiques et les profils des horizons 
g^ologiques avec un contraste am^ior^ par rapport k Tart anterieur. 

Pour une autre variante de I'invention k I'dt^e 105, on choisit comme indice' 
de courbure de Thorizon au point de coordonn^ spatiales i, j et de coordonn^e 
temporelle ou de profondeur k une fonction croissante des rdsidus comportant les 
coefficients ayjcet bij.k ddterminds h I'dtape 103. 

On pent notamment utiliser comme indicateur ou indice de courbure le critfere 
de minimisation C«(id,k) = X(Pu.pq.k)" ^^s6 pour calculer les coefficients 
param^triques aijjcetbyjc d6termin& k Mape 103. 

On peut egalement utiliser sans sortir du cadre de la pr€sente invention, toute 
autre fonction croissante des rdsidus sur le voisinage conditionnel Njj.k(r). 

En effet, dans le cas oii I'horizon est parfaitement plan au voisinage du point 
de coordonndes spatiales i j et de coordonnde temporelle ou de profondeur k. toute 
fonction croissante des r&idus sera ^gale k z6ro et I'indice de courbure locale sera 
par consequent Egalement dgal k zj6to. 



Au contraire, dans le cas ot Thorizon au point de coordonn^es spatiales i,j et 
de coordonnee temporelle ou de profondeur k n'est pas plan, rindice de courbure 
calculi k I'dtape 105 du pTOc6d6 selon Tinvention fourait une indication ou une 
mesure de la difference entre rhorizon consider^ et le plan localement osculateur au 
point de coordonn^es spatiales i J et de coordonnee temporelle ou de profondeur k. 

^invention permet dgalement par le calcul de Tindice de courbure et sa 
representation graphique sur iin ecran d'amdliorer le contraste des mesures sismiques 
fournies par les geophones sous forme de matrice sismique tridimensionnelle et 
d*6viter ainsi des mesures sismiques additionnelles jugees au voisinage de 
discontinuitds ou de failles gdologiques. 

La mise en oeuvre pr^ferfe d'un procdde selon Tinvention consiste h utiliser 
un dispositif pour la mise en oeuvre du proc€de contenant un produit-programme 
d'ordinateur ex&utant les etapes du proc^de selon Tinvention ddcrites en reference 
aux figures 1^6. 

A cet effet, le dispositif selon llnvention est muni de moyens de 
memorisation appropries pour m^moriser et calculer les elements successivement 
utilises dans le procede selon Tinvention. 

De preference, la mise en oeuvre du procede selon Tinvendon s'effectue sur 
une station de travail de calculs sismiques en chargeant un produit-progranmie selon 
rinvention comportant des elements de codes de programme pour executer les etapes 
du procede selon Tinvention. 

Le produit-programme peut comporter des elements de codes de progranmie 
relatifs aux etapes 100, 101, 102, 103 et 104 decrites en reference k la figure 6. 

Alternativement, le produit-prograname selon Tinvention peut comporter les 
elements de codes de programme pour executer les etapes 100, 101, 102, 103 et 105 
decrites en reference a la figure 6. 

Selon d'autres variantes, un produit-progranmie d'ordinateur selon I'invention 
peut egalement comporter des elements de codes de prograname pour executer au 
moins Tensemble des etapes de procede decrites en reference k la figure 6, 

Uinvention decrite en reference k plusiexirs modes de realisation n'y est 
nuUement limitee, mais couvre au conUaire toute modification de forme et toute 




variante de mise en oeuvre, incluant ou non une presentation de r€sultats 
interm^diaires sur €cran de visualisation* 



REVENDICATIONS 

1. Pvoc6d6 pour determiner localement la forme des horizons gdologiques, 
,dans— LequeL^oa. constniit une . foncJiQa..,xmtimie^^ijjcQj^ 

approximation des traces sismiques discrfetes d*une matrice sismique 
5 tridimensionnelle, ladite fonction SyjcCt) ^tant d6sign6e comme "trace sismique locale 
continue", comportant les Stapes suivantes : 

a) , utiliser comme decalage optimal de deux traces sismiques locales 
continues voisines, la valeur de decalage rendant maximale leur fonction de 
correlation, ce ddcalage optimal n'dtant pas contraint k 6tre un multiple entier du pas 

10 d'dchantillonnage vertical ; 

b) . retenir comme voisinage conditionnel d*une trace sisnoique locale continue 
"centrale" Sij,k(t) le sous-voisinage consistant en des traces voisines Spq,k(t) 
correspondant h, des ddcalages optimaux hij,pq,k associes k des correlations 

15 Rij.pq.k(hii.pq,k) superfeures k un seuil predetermine compris entre 0 et 1 ; 

c) . definir pour chaque trace sismique locale continue Spq,k(t) du voisinage 
conditionnel une valeur de residu par rapport k ladite trace sismique locale continue 
"centrale" Sij,k(t) comportant des coefficients parametriques ; 

d) . determiner les coefficients parametriques ay^k et by^ par nainimisation 
20 d'ensemble des residus sur le voisinage conditionnel. 

2. Procede selon la revendication 1, dans lequel les valeurs de residu de 
retape c) sont des valeurs absolues de forme differentielle parametrique, par exemple 
I ay,k .(p-i) + bij,k .(q-j) - hij.pq.k I , ou ay^k et byjc sont les parametres et hy,pq.k est le 
decalage optimal retenu entre la trace sismique locale continue Sy,k(t) "centrale" et 

25 une trace locale sismique locale continue Spq,k(t) voisine. 

3. Procede selon la revendication 1 ou la revendication 2, dans lequel la 
minimisation d'ensemble des residus Pij,pq,k sur le voisinage conditionnel comporte 
une minimisation d'une somme de puissances des residus, de la forme 



C ^ (io,k) = X(Aj.pq.k)" » ^®st une puissance superieure k 1. 
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4. Pioc6d6 selon la revendication 1, dans lequel les coefficients aij.ket bu.k d^ermin^s 
k l'6tape b) sont utilises pour d€finir un vecteur unitaite NCiji) de coordonn€es 
N''(ij,k) , N5'(i j^c) , N*(ij,k) , par exemple de la forme: 



et dans lequel on choisit ce vecteur unitaire conrnie vecteur unltaire normal 
(orthogonal) k I'horizon du nceud (i, j, k). 

5. Proc^d6 selon la revendication 4, dans lequel on utilise I'ensemble des 
vecteurs unitaires normaux N(i. j, k) pour d^finir le champ des vecteurs unitaires 
norraaux, en vue de I'affichage sur 6cran des profils des horizons sisraiques. 

6. Proc^6 selon I'une quelconque des revendications 1^3, dans lequd on 
choisit comme indice de courbure une fonction croissante des r&idus comportanit les 
coefficients paramdtriques d^termln& k rstape d). - 

7. Dispositif pour la mise en ceuvre du precede selon Tune quelconqu6 des 
revendications 1 k 4, comportant des moyens pour utiliser comme decalage optimal 
de deux traces sismiques locales continues voisines la valeur de ddcalage rendant 
maximale leur fonction de coixdlation, des moyens poiir retenir comme voisinage 
conditionnel d'une trace sismique locale continue "centrale" SyjcCt) le sous-voisinage 
consistant en des traces voisines Spqjc(t) conespondant h des d&alages optimaux 
associes k des correlations Rij.pqjc(h) sup^rieures k un seuil prdddtermind compris 
entre 0 et 1, des moyens pour d^finir pour chaque trace sismique locale continue 
Spqjc(t) du voisinage conditionnel une valeur de rdsidu par rapport k ladite trace 
sismique locale continue "centrale" Sijjc(t) comportant des coefficients param^triques 
et des moyens pour d^erminer les coefficients paramdtriques par minimisation 
d'ensemble des r^sidus sur le voisinage conditionnel. 

8. Dispositif selon la revendication 7, comportant des moyens de 
memorisation et des moyens de visualisation de paramfetres sismiques determines k 
I'aide du procede selon I'lme quelconque des revendications Ik 6. 



N"(iJ,k) = 




Ny(ij,k) = 



■v/Kk)'+(b«.k)'+l 
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9. Produit-programme d^ordinateur, compoitant des 6l6ments de code de 
programme pour executer les dtapes du proc€d6 selon Tune quelconque des 
revendications 1 a 5, lorsque ledit programme est execute par un ordinateur. 

10. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour executer les Stapes du proc€d6 selon la revendication 6, lorsque 
ledit programme est tx6cut6 par un ordinateur. 
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